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Abstract 

We present an ab initio pseudopotential calculation within the local density aproxima- 
tion of the electronic structure and the full set of elastic constants of GaAs and ZnSe.We 
compare with the results established experimentally, to determinate the importance of 
the core corrections in Zn and Ga and we find that we must take core corrections inlo 
account to obtain accurate results. 

Resumen 

En este trabajo presentamos una determinación de la estructura de bandas y de las cons¬ 
tantes elásticas del ZnSe y el GaAs utilizando pseudopotenciales ab initio y conservadores 
de la norma dentro del formalismo del funcional de la densidad. Lo^ resultados obtenidos 
se comparan con los valores obtenidos experimentalmente para determinar la importancia 
de las correcciones de ”core” , en el Zn y el Ga. Nuestros resultados confirman que las 
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correcciones deben tenerse en cuenta cuando tratamos con compuestos formados a partir 
de Zn y Ga. 

I Introducción. 

En la década pasada comenzaron a desarrollarse los denominados cálculos ab initio (cálculos 
desde primeros principios), en los que no se incluye ningún tipo de ajuste a partir de datos 
experimentales. Entre las propiedades que se pueden obtener de esta forma se encuentran: las 
estructuras cristalinas, el espectro de fonones, las constantes de red, los módulos de compresión 
y cizalla y otras propiedades estáticas y dinámicas. Este desarrollo ha abierto muchas posibili¬ 
dades para estudiar la materia condensada, puesto que en la actualidad es incluso posible una 
buena predicción de las propiedades que no pueden obtenerse de forma experimental debido a 
las características del sistema [2]. 

En este trabajo determinaremos , partiendo de primeros principios, los módulos de elastici¬ 
dad y las estructuras de bandas de dos semiconductores del tipo III V y II-VI que cristalizan 
en la estructura zinc-blenda. En concreto nos centraremos en el GaAs , cuyas propiedades 
ya han sido ampliamente estudiadas tanto desde el punto de vista experimental como desde el 
punto de vista teórico, y en el ZnSc no tan estudiado como el compuesto anterior. 

Utilizamos la teoría del funcional de la densidad en la aproximación de la densidad local 
[1], LDA (Local Density Approximation). Consideraremos además en nuestro planteamiento 
inicial, y a la hora de realizar los cálculos, la inclusión de las denominadas correcciones no 
lineales de core . Estas correcciones consisten en un adecuado tratamiento en la energía de 
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canje y correlación, de los electrones de la capa 3 d (electrones del core) que se encuentran r íy 
próximos energéticamente a los electrones de la capa de valencia. Para el ZnSe, los cálculos 
electrónicos efectuados muestran que las bandas M del Ztl están sólo a 7 eV aproximadamente 
por debajo del máximo de la banda de valencia presentando un pequeño grado de dispersión [8], 
mientras que para el GaAs las bandas 3 d del Ga están aproximadamente 2.7 eV por debajo de 
dicho máximo y no presentan casi dispersión [4], Por tanto, para los semiconductores formados 
a partir del Ztl el error en el que se incurre tratando los electrones 3 d como estados de core 
congelados (frozen core States), es aún más importante que para los semiconductores formados 
a partir de Ga [10]. En consecuencia, para tener en cuenta los efectos de los electrones 3 d del Zn 
y del Ga , es necesario a su vez la introducción de las denominadas correcciones de correlación 
e intercambio no lineales , (dada la forma no lineal de la energía de canje y correlación), para 
obtener una descripción razonable de las propiedades estructurales para los dos semiconductores 
estudiados [11]. 

El método de cálculo utilizado para determinar las constantes elásticas fue desarrollado por 
Nielsen y Martin [12] cuando estudiaron las constantes elásticas para el Ge, el Si y el GaAs 
usando los cálculos de la tensión macroscópica y de la fuerza sobre los átomos en el sólido. 
Basándonos en cálculos previos efectuados para materiales similares [2,3], podemos concluir 
que nuestros resultados son igualmente exactos y pueden ser usados para obtener la estructura 
de bandas electrónica y las constantes elásticas tanto del GaAs como del ZnSe. 

En la sección II describiremos de forma breve el método de cálculo. En la sección III 
mostraremos los resultados obtenidos para las propiedades electrónicas y las constantes elásticas 
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del GaAs y el ZnSe y los compararemos con los de otros trabajos anteriores. Finalmente en 
la sección IV presentamos las conclusiones obtenidas. 

II Método de cálculo 

En este estudio hacemos uso de la teoría del funcional de la densidad dentro de la aproxi¬ 
mación de la densidad local , LDA, [1] empleando un desarrollo en ondas planas y pseudopoten- 
ciales no locales y conservadores de la norma del tipo de Bachelet, Hamann y Schlüter [4], junto 
con la forma de Ceperly y Alder para el canje y la correlación [5] tal y como fue parametrizada 
por Perdew y Zunger [6]. Se tendrán en cuenta además las denominadas correcciones no 
lineales de core, siguiendo el formalismo propuesto por Louie [7]. 

Las propiedades elásticas se determinan calculando las componentes del tensor de tensiones 
para pequeñas deformaciones usando el método desarrollado por Nielsen y Martin [12]. Para 
una deformación pequeña en la direccipn (100) las constantes elásticas C\\ y Cu se obtienen 
de las relaciones en el límite armónico C\\—(T\lci y Cu—cj^I^x-í donde cr t y e,- representan 
respectivamente, la tensión y la deformación aplicada usando la notación de Voigt [14]. La 
tensión macroscópica en el sólido se calcula para una pequeña deformación usando el Teorema 
de la tensión [12] y las fuerzas sobre los átomos se obtienen a partir del Teorema de Hellmann- 
Feynman [16,17]. Para una deformación e 4 (deformación uniaxial en la dirección [111]) existen 
desplazamientos internos de las subredes y las posiciones atómicas en la celda unidad no quedan 
determinadas únicamente por simetría. Para estos casos Kleinman [18] definió un parámetro 
de deformación interna £ que describe el desplazamiento de los átomos. Nielsen y Martin [12] 
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formulan dos cálculos independientes de la tensión y la fuerza que determinan tres magnitudes 
independientes C44, £ y la frecuencia del fonón óptico R Con ello se obtiene la siguiente 
relación tensión-deformación: 

= [C° 4 - n-'<S>(ía 0 /4) 2 )t 4 = C 4A e< 

donde C 44 es la constante elástica en ausencia de desplazamientos internos y 4> es la constante 
de fuerza, que es igual a puf*, donde p es la masa reducida, a 0 es la constante de red en el 
equilibrio y u es el volumen de la celda unidad. 

III Resultados 

A. Constantes Elásticas 

En orden a obtener resultados bien convergidos las funciones de onda se expanden en 
términos de ondas planas con una energía cinética de corte de 24 Rv para ambos semicon¬ 
ductores [13]. En la Fig.l mostramos la presión isotrópica del ZnSe en función de la energía 
de corte. En dicha figura puede observarse como a partir de 24 Ry la presión isotróprica se 
mantiene constante, lo cual indica que esta energía de corte es suficiente para tener resulta¬ 
dos bien convergidos [12]. La integración sobre la zona de Brillouin se realiza con suficiente 
aproximación utilizando un conjunto de diez puntos especiales [15] . Las constantes de red 
del equilibrio se determinan ajustando la energía total a la ecuación empírica de estado de 
Murnaghan [19]. La energía total en función del volumen por átomo para el caso del ZnSe se 
muestra en la Fig.2. 
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Para obtener las constantes elásticas y la frecuencia del fonón para el ZnSe usamos una 
deformación de C\=± 0.002 en la dirección (100) para obtener C\\ y Ci 2 . Para determinar ¿jp, 
C44 y £ aplicamos un pequeño desplazamiento u = ±0.002 A a lo largo de la dirección (111) 
y una deformación e 4 =± 0.002 en la misma dirección. 

Nuestros resultados obtenido utilizando las correciones de core para las constantes elásticas 
en el caso del ZnSe son Cn=1.125 Mbar, Ci 2 =0.501 Mbar y 644 = 0.61 Mbar 

El módulo de compresión está definido como B = (C n + 2C i2 )/3 y la constante de cizalla 
efectiva como C s = (Cu — C i2 /2. 

Para el GaAs la deformación aplicada a lo largo de la dirección (100) para obtener C\\ 
y Ci 2 es €i=± 0.004. Mientras que para determinar cjp, C 44 y ( se aplica un desplazamiento 
« = ±0.004 Á y una deformación e 4 =± 0.004 a lo largo de la dirección (111). Los resultados 
obtenidos considerando y sin considerar correcciones de core en el Ga se muestran en la Tabla 
I y para el ZnSe en la Tabla II. 

De la observación de las tablas de resultados podemos extraer que los valores obtenidos para 
el módulo de compresión y la frecuencia up concuerdan notablemente con los valores experi¬ 
mentales [23,24] y [27,28]. Las constantes elásticas calculadas comparan bien con los resultados 
obtenidos mediante scattering de Brillouin [22,26], esta concordancia con los resultados expe¬ 
rimentales no existe si se desprecian las correcciones de core. 

B. Estructuras de Bandas 

Continuando con el estudio de las propiedades de volumen de los dos semiconductores 
tratados, se procedió al cálculo de sus respectivas estructuras de bandas. Mostrándose los 
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resultados obtenidos para el GaAs y el ZnSe con correcciones de core en el Ga y el Zn en las 

Figs.3 y 4, . 

Teniendo en cuenta que el método utilizado en nuestro estudio no describe de forma exacta la 
energía de los estados excitados aunque nos proporciona una descripción cualitativa de éstos, de 
la observación de las estructuras de bandas obtenidas para el GaAs y para el ZnSe, incluyendo 
correcciones de core en el Ga y en el Zn, vemos como el grado de concordancia entre los 
resultados anteriores conocidos [8,20,30] y los obtenidos en nuestro trabajo es excelente para 
las bandas de valencia asi como para la forma de las bandas de conducción (Es conocido que la 
teoría del funcional de la densidad solo proprociona información, estrictamente correcta, sobre 
el estado fundamental, subestimando generalmente el gap). 

IV Conclusiones 

En este trabajo y usando métodos desde primeros principios hemos obtenido energías totales 
con un elevado grado de convergencia, fuerzas y tensiones en el GaAs y el ZnSe para llegar 
a los módulos de elasticidad, las frecuencias del fonón TO( F) y las estructuras de bandas de 
estos materiales. Para los cálculos anteriores hemos usado un método autoconsistente con 
pseudopotenciales ab initio en el formalismo LDA con una base en ondas planas. En base a 
la observación de los resultados experimentales ya tabulados y los de otros cálculos en los que 
no se tenían en cuenta las correcciones de core, se llegó a la conclusión de que los electrones 
de la capa 3 d tanto del Ga como del Zn no habían sido tratados de forma adecuada. Por lo 
tanto, en nuestro estudio usamos las denominadas correcciones de correlación e intercambio no 
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lineales que permiten un adecuado tratamiento de estos electrones. A la vista de los resultados 
obtenidos concluimos que estas correcciones de core permiten una mejor descripción de las 
propiedades estructurales de los dos compuestos, GaAs y ZnSe , estudiados. 
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Tablas 


TABLA I. (a) Constante de red a 0 , constantes elásticas C XJ del GaAs con correcciones de 
core en el Ga, la constante elástica CJ 4 , la frecuencia del fonón transversal óptico en I\ ir , 
el módulo de compresión B. la constante de cizalla efectiva C s y el parámetro de deformación 
interna . (b) El mismo estudio sin correcciones de core. 


Ref. 

a 0 

(A) 

C\\ 

(Mbar) 

C\2 

(Mbar) 

(Mbar) 

C°„ 

(Mbar) 

*T 

(THz) 

B 

(Mbar) 

c, 

(Mbar) 

í 

Expt. [ 21 ] 

5.653 









Expt. [ 22 ] 


1.126 

0.571 

0.60 






Expt. [23] 






8.013 




Expt. [24] 







0.769 



Referencia [3] 


1.38 

0.55 

0.66 






Referencia [12] 

5.55 

1.23 

0.53 

0.62 


8.09 

0.73 


0.4S 

Este trabajo (a). 

5.643 

1.125 

0.501 

0.61 

0.761 

8.04 

0.709 

0.312 

0.314 

Este trabajo(b). 

5.387 

1.533 

0.685 

0.82 

1.036 

8.9 

0.96 

0.424 

0.325 








0.995 1 




1 Resultado obtenido ajustando mediante la ecuación empírica de Murnagham. 
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TABLA II.(a) Constante de red a 0 , constantes elásticas C tJ del ZnSe con correcciones de 
core en el Zn, la constante elástica Cj 4 , la frecuencia del fonón transversal óptico en T, i/r, 
el módulo de compresión B. la constante de cizalla efectiva C s y el parámetro de deformación 
interna £. (b) El mismo estudio sin correcciones de core. 

Ref. üq C\\ C \2 C 44 C44 vy B C s £ 

(Á) (Mbar) (Mbar) (Mbar) (Mbar) (THz) (Mbar) (Mbar) 

Expt. [25] 5.667 

Expt. [26] 

Expt. [27] 

Expt. [28] 

Referencia [ 8 ] 5.409 

Referencia [29] 5.53 

Este Trabajo(a) 5.632 
Este trabajo(b) 5.183 


0.872 0.524 0.392 


0.873 

0.514 

0.43 

0.56 

1.4 

0.79 

0.758 

1.009 


0.667 

0.283 


6.2 

0.635 

0.175 

0.377 

7.9 

1.01 

0.32 

0.37 


0.872 2 




2 Resultado obtenido ajustando mediante la ecuación empírica de Murnagham. 






Presión isotrópica del ZnSe en función de la energía de corte para una constante de red 
a 0 = 5.565 A. 



VOLUMEN/ATOMO (A) 


Fig.2 

Energía total del ZnSe en función del volumen por átomo a una energía de corte de 24 Ry. 
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Estructura de bandas del GaAs calculada con una energía de corte de 24 Ry y una constante 
de red a Q = 5.643 A. 



Estructura de bandas del ZnSe calculada con una energía de corte de 24 Ry y una constante 


de red a 0 = 5.632 A. 
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